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Abstract 

Electron magnetohydrodynamic (EMHD) turbulence in two dimensions is studied via high- 
resolution numerical simulations with a normal diffusivity. The resulting energy spectra asymp- 
totically approach a law with increasing Rb, the ratio of the nonlinear to linear timescales 
in the governing equation. No evidence is found of a dissipative cutoff, consistent with non-local 
spectral energy transfer. Dissipative cutoffs found in previous studies are explained as artificial 
effects of hyperdiffusivity. Relatively stationary structures are found to develop in time, rather 
than the variability found in ordinary or MHD turbulence. Further, EMHD turbulence displays 
scale-dependent anisotropy with reduced energy transfer in the direction parallel to the uniform 
background field, consistent with previous studies. Finally, the governing equation is found to yield 
an inverse cascade, at least partially transferring magnetic energy from small to large scales. 
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I. INTRODUCTION 



Turbulence plays a crucial role in a wide variety of geophysical and astrophysical fluid 
flows. In this paper we present results on a particular type of plasma turbulence in which 
the flow is entirely due to the electrons, with the ions forming a static background. The 
equation governing the electrons' self-induced magnetic field is then 

^ = - V X [J X B] + R-^' V^B, (1) 

where J = V x B, and Rb = oB^jnec^ with a the conductivity, a measure of the 
field strength, n the electron number density, e the electron charge, and c the speed of 
light. See for example who derived this equation in the context of magnetic fields in the 
crusts of neutron stars. More generally though, it is applicable in many weakly collisional, 
strongly magnetic plasmas, so other applications could include the Sun's corona or the 
Earth's magnetosphere. 

Turbulence governed by ([1]) is known as electron MHD (EMHD), Hall MHD, or whistler 
turbulence. Based on its (at least superficial) similarity to the vorticity equation governing 
ordinary, nonmagnetic turbulence, 

^ = V X [u X w] + Re-^ VV, (2) 

where now w = V x u, [l] argued that ([1]) would initiate a turbulent cascade to small 
lengthscales, thereby accelerating neutron stars' magnetic field decay beyond what ohmic 
decay acting on large lengthscales could achieve. They suggested in particular that the 
turbulent spectrum would scale as A;~^, with a dissipative cutoff occurring at A; ~ Rb- 

However, there are also some quite fundamental differences between equations ([1]) and 
([2]). In ([2]) the dissipative term contains more derivatives than the nonlinear term, so on 
sufficiently short lengthscales the dissipative term will always dominate. In contrast, in ([1]) 
the two terms both contain two derivatives, so it is conceivable that the nonlinear term will 
always dominate, even on arbitrarily short lengthscales. As pointed out by [2|, one obtains 
a dissipative cutoff only if one assumes that the cascade is local in Fourier space, coupling 
wavenumbers only to their immediate neighbors. 

Indeed, whether the coupling is local or not is another important difference between ([1]) 
and ([2]). In ([2]) it is at least predominantly local; for example, very small scale structures see 
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the largest structures as an essentially uniform background flow that simply advects them 
along, but without altering the nature of the small scale turbulence. In contrast, in ([1]) there 
is no such translational invariance; adding even an exactly uniform background field alters 
the dynamics of the small scale structures (as we will show in detail below). Similarly, in 
classical MHD turbulence, adding a background flow has no essential effect, but adding a 
background field does. 

In this paper we present high-resolution numerical simulations of ([1]) in a two-dimensional 
periodic box geometry, designed specifically to address such questions as whether there is a 
dissipative cutoff or not, and whether the coupling is local or not. In contrast to previous 
simulations j^, ^, 0, [g], Q, [s], we do not employ hyperdiffusivity, which would of course 
disrupt this feature that the two terms in ([T]) have the same number of derivatives, and 
hence introduce an artificial dissipative cutoff. We also consider the question of whether 
([1]) is capable of yielding an inverse cascade, and find that magnetic energy can be at least 
partially transferred from small to large scales. 



II. EQUATIONS 

For two-dimensional fields, we may decompose B as 

B = Bp + = V X (ae^) + be^. (3) 
where a and b depend only on x, y, t, but not z. Equation ([1]) then yields 

— = R]^^V^a - {ttyb^ - a^by), (4) 



and 



db 

— = R^^V^b + {ayV^a, - a^V^ay), (5) 



where subscripts indicate derivatives. 

Continuing our comparison of equations ([1]) and ([2]), it is instructive to note also that in 
([2]) we would only have u = Up = V x (^Pe^), yielding 

= i?e-iV% - [v^, (V^vl/), - (V^vl/),]. (6) 

Any additional ut = vez would simply be advected by \E' as a passive scalar, but without 
any influence back on \I'. This difference between (jlj) and ([5]) on the one hand, and ([6]) on 
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the other merely reflects once again some of the differences between ([T]) and ([2]), in this case 
the lack of translational invariance in ([1]). 

We solve (jl]) and by expanding a and b in Fourier series in x and y, and using stan- 
dard pseudospectral techniques for the evaluation of the nonlinear terms, with dealiasing 
according to the 2/3 rule. The code employs the FFTW library [9] to achieve massive par- 
allelisation on a suitable supercomputer. The time integration is done using a second order 
Runge-Kutta method. We performed a variety of runs, typically employing 64 processors, 
with the highest extending to A; = 682 in Fourier space, corresponding to = 2048 colloca- 
tion points in real configuration space. Because of the two derivatives in the nonlinear term, 
the required timesteps are unfortunately very small, roughly proportional to 1/ (A^^). Values 
as small as ~ 3 x 10~^ were used, requiring O(IO^) timesteps in total to reach t = 0.2. 

A. Initial Conditions 

Since we are interested in freely decaying rather than forced turbulence, we need to 
carefully consider the nature of our chosen initial conditions. We will present results for 
three different sets of runs. 

First, to study homogeneous forward cascades, we start off with random 0(1) energies in 
all Fourier modes up to = ^/fc^ + ky = 5, making sure that the poloidal a and toroidal b 
components have comparable amounts of energy. After initialisation the overall amplitude 
of the field is rescaled to ensure that the rms value of |B| = 1 at t = 0. 

Second, to study nonhomogeneous forward cascades, we start off with the same initiali- 
sation as above, but now add a uniform field Ce^, where C = 1, 2 or 4. Note though that 
such a uniform field cannot be represented by an expansion of the form ([3]), at least not if 
a and b are to be periodic in x and y. Instead, this field is simply added in to ([3]) directly, 
resulting in suitably modified equations (jl]) and ([5]). 

Third, to explore the possibility of inverse cascades, we return to the C = case without 
any large scale magnetic field and now inject energy into modes in the range 10 < k < 20. 
The question then is how much of this initial energy moves to k < 10, and how much moves 
to A; > 20. 

Finally, for all three sets of results, each individual run was repeated with a number 
of different random initial conditions, to ensure that the results presented here are indeed 
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representative. 



B. Ideal Invariants 



Equations (jlj) and ([5]) also have some useful associated diagnostics, corresponding to 
quantities that are conserved in the ideal, i?^^ 0, limit. Specifically, we have equations 
for the energy and the magnetic helicity. 



dt2 



-R 



B 



d 1 
di2 



A-BdV = -R^^ 



B • J d\/. 



(7) 



(8) 



where A is the vector potential, defined by B = V x A. Note though that in the presence 
of a uniform background field, helicity is not even defined lOj, let alone conserved. 

These two equations are valid in both 2 and 3 dimensions. In 2 dimensions only, we have 
the additional quantity of the mean squared magnetic potential, known as anastrophy 



d 1 
dt2 



a2 dV = -R^^ 



ai + al) dF, 



(9) 



which is in some ways perhaps analogous to enstrophy, which is also defined only for ([6]) 
in 2D, but not for ([2]) in 3D. However, anastrophy is not the same as enstrophy, and there 
does not appear to be any reason why conservation of anastrophy would necessarily imply 
an inverse cascade in the way that conservation of enstrophy forces inverse cascades to exist 
in 2D hydrodynamic turbulence. 

In addition to the physical insight that they yield into the nature of the Hall nonlinearity, 
these various integrated quantities also offer useful diagnostic checks of the code. Reassur- 
ingly, we found that all of them (except helicity in a uniform field of course) were satisfied 
to within 0.1% or better by all of our runs. 

III. RESULTS 



A. Large-scale initial conditions 



The energy spectrum, E^, o: 
the familiar Kolmogorov law 



' 3D hydrodynamic turbulence is characterised by oc k 

In 2D, the conservation of enstrophy forces an inverse 
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cascade which leads to a much steeper spectrum (x k " with a > 3. In the case of 2D 
MHD turbulence, the spectral energy transfer rate is reduced which leads to a flatter energy 



3- 



Recent studies of EMHD 



spectrum oc the Iroshnikov-Kraichnan spectrum 

turbulence have found, via methods which all employ hyperdiffusivity, a 5/3 Kolmogorov 
spectrum for small scales k(L> 1, equivalent to A; > 0{Rb), and a steeper 7/3 spectrum for 
longer wavelengths j^, 4, 5, 6, 7|. 

In the upper plot of Figure [H we show the poloidal and toroidal energy spectra of our 
solutions for Rb = 10, 30, 100, 300, 1000 & 3000, evolved to a time t = 0.2. The energy 
spectra have been stationary since approximately t = 0.12 and time averaging between 
0.12 and 0.2 reveals an identical spectrum and no further information. We interpret this 
to mean our simulations are resolved and evolved to a suitable time for inspection of the 
quasi- stationary cascade. 

Both poloidal and toroidal energy spectra start out much the same at low k and then 
lower Rb spectra smoothly drop off with increasing k whilst higher Rb spectra maintain 
a linear gradient in the log-log plot. Transfer of energy to higher k is then more efficient 
at higher k. The spectra are asymptotically approaching an energy spectrum Ek oc k~'^, 
where u = 2.5 ± 0.1. In the lower plot of Figure [1], we show compensated energy spectra to 
show this approach to k~^/'^Ek = 1 with increasing Rb- Our value of u is not compatible 
with —2 predicted by |l| for 3D EMHD turbulence. Their prediction was calculated using a 
phenomenology based on Kraichnan's arguments (the whistler effect) which has been shown 
by to have little effect on the energy spectrum of 2D EMHD, rather it is thought to 
influence the subtle properties of the cascade like anisotropy which we discuss further in the 
next section. The spectral index we find is much more compatible with u = 7/3 found by 
jl] via numerical simulation for the kde > 1 regime. 

None of the spectra show any sign of a dissipative cutoff. By definition, the dissipation 
scale should occur when the local value of Rb is 0(1) in equation [H It is unclear though 
when this occurs since the definition oi Rb does not involve length scales. If the coupling 
is purely local in wavenumber, then this definition does involve length scales after all, since 

the Bq that should be used is the field at that wavenumber only, rather than the total field. 

i. 

That is, according to the definition of where this argument was first developed, we have 

R's = Rb{B'/B) (10) 
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where the primed quantities are the small-scale local values and the unprimed the large-scale 
global. If we now suppose a energy spectrum, then B' / B ~ and so R'^ is reduced 
to 0(1) when k ~ R^^^ ■ So, at -R^ = 1000 we would see a dissipative cutoff at A; ~ 250. 
We do not see a dissipative cutoff at this scale. We can reconcile this by realising that this 
argument crucially depends on the coupling being local in Fourier space: if this does not 
hold then R'^ = Rb and there is simply no definite dissipation scale, as we find here. 

We conclude therefore that the nonlinear term is able to dominate at all length scales 
and the coupling is non-local in Fourier space. This is not entirely unexpected, since both 
the terms in ([T]) contain two derivatives, but it is in contrast to previous studies. It should 
be emphasized that other authors have been unable to properly address the question of a 
dissipative cutoff since hyperdiffusivity has masked the effect of the nonlinear term at high 
k. This may also explain the discrepancy between the spectral index of 5/2 we find and 
previous values of 7/3. 

In Fourier space then, EMHD turbulence bears a strong resemblance to ordinary MHD 
turbulence. We would like to know if this resemblance carries over into real configuration 
space. In Figure [2], we show the three component fields at three times; the initial fields at 
t = (top row), intermediate fields at t = 0.1 (middle row) and fully developed turbulent 
fields at t = 0.2 (bottom row). In classical and MHD turbulence, a fully developed turbulent 
field in real configuration space would bear no resemblance to the initial field, but here the 
fields are much more structured and fully developed turbulent fields strongly resemble initial 
fields. This appears to be a unique characteristic of decaying EMHD turbulence. 

We would also like to address the energy decay of the field, with particular respect to any 
dependency of the decay rate on the value of Rb- Ref. [3| in the first study of 2D EMHD 
reported that the energy dissipation rate is independent of the value of the dissipation 
coefficient, represented hj Rb here. In contrast to this, we find that the energy decay is 
much slower at higher Rb, as plotted in Figure O 



B. Large-scale initial conditions in the presence of a background field 

EMHD turbulence, like classical and MHD turbulence, is isotropic when allowed to freely 
decay. In the presence of a background flow, classical turbulence remains isotropic since 
it is locally coupled in Fourier space. Small-scale structures are simply advected along by 
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the large-scale flow, whether or not that has a uniform background contribution. Numerical 
simulations of MHD turbulence have found it to be strongly anisotropic in the presence 



of a background field 
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14|. This has been attributed to the excitation of Alfven waves 



which preferentially propagate parallel to the external magnetic field and hinder the cascade 
process perpendicular to the external field. 

In EMHD turbulence, recent numerical studies employing hyperdiffusivity have 
revealed similar strongly anisotropic behaviour. This can only be the result of asymmetry in 
the nonlinear spectral transfer process relative to the external magnetic field. In the context 
of local energy coupling in Fourier space, mediation by whistler waves has been proposed 
as the only way this asymmetry could be achieved a] with the method detailed in ref. 15|. 
The spectrum of 2D anisotropic EMHD turbulence has also been shown to exhibit a linear 
relationship with an external magnetic field Q]. 

In order to understand how hyperdiffusivity has affected previous studies we have intro- 
duced a background field into the governing equations as discussed above and calculated 
solutions for Rb = 100 & 300 at a spatial resolution of 512^ points in real configuration 
space. We present our results in Figure |H Across the top row, we show 2D energy spectra 
for Rb = 100 with C = 0, 1, 2 & 4. In the isotropic case with no background field, i.e. 
C = 0, energy is evenly distributed between x and y, as indicated by circular contours. In 
the case of C = 1 we find energy transfer to larger k has been suppressed in the x direction, 
parallel to the background field. EMHD turbulence has become anisotropic in the presence 
of a uniform background field with normal diffusivity. The effect becomes more pronounced 
for C = 2 and C = 4. The evolution of modes parallel and perpendicular to the field is 
clearly different, as the spectral cascade in the parallel wavenumbers is clearly suppressed. 
This suppression has been attributed to excitation of whistler waves, which act to weaken 
spectral transfer along the direction of propagation jsjl . Across the bottom row, we show the 
corresponding fields. For C = 0, the field is isotropic, but as the value of C is increased, 
structures are stretched in the x direction corresponding to increasingly inhibited energy 
transfer in x but not in y, perpendicular to the field direction. 

This result is in agreement with [7] who found scale-dependent anisotropy in numerical 
studies of 3D EMHD turbulence employing hyperdiffusivity. Our simulations also support 
the linear relationship between 2D EMHD turbulence and strength of external magnetic 
field found by j^. At Rb = 300 we reassuringly find the field is more anisotropic. 
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It is worth noting here that since ([T]) is scale invariant - we can apply the equation over 
the whole of a system, or just a small section, with Rb unchanged - if you take a very small 
box, then this box will see the large-scale field as a background field, and therefore one would 
expect the smallest scales in the system, for example a neutron star, to be anisotropic. 

C. Intermediate-scale initial conditions 

In classical turbulence, the exchange of energy and enstrophy Q is coupled in Fourier 
space according to 

dt~ " dt' ^^^^ 

hence energy injected at intermediate scales experiences a transfer to both higher and lower 
wavenumbers in order to satisfy this coupling and simultaneously conserve energy and en- 
strophy. This is the inverse cascade of energy to lower k (larger scales) [l6|. In MHD 
turbulence, energy and magnetic helicity are coupled in the same way and an inverse cas- 
cade occurs in order to simultaneously conserve these two quadratic ideal invariants. To 
investigate whether an inverse cascade occurs in decaying EMHD turbulence, we inject en- 
ergy over the wavenumber range 10 < A; < 20 as detailed above and evolve the magnetic 
field. 

In Figure [5] we show the solution for Rb = 1000 at various times. To fully resolve the 
solution in a reasonable amount of computational time we have chosen to evolve the field 
to t = 0.2 at a resolution of 2048^ real space points, then to t = 1.0 with 1024^ points 
and finally to t = 15.0 with 512^ points. For this reason the energy spectra have different 
extents in Fourier space at the different times t = 0.0 (dashed line), 0.1, 1.0, 3.0, 6.0, & 
15.0. We include the full information, rather than just cutting off the plot at k = 100, 
to demonstrate that our solutions are indeed fully resolved. The spectra show that energy 
is clearly transferred to < 10 in an inverse cascade, with the spectral peak shifting to 
k = 2 but not maintaining the same amplitude. Some energy has also transferred to > 20 
resulting in an overall spectrum comparable to a forward homogeneous cascade at late time. 
At t = 0.1 the spectrum has a spectral index of -5/2 for k > 10. At late time this has 
steepened to ~ —3. Between t = 12 and t = 15, the spectral index has stabilised. The 
inverse cascade phenomenon becomes less pronounced at lower values of Rb with no inverse 
cascade at all below Rb = 300. 
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In Figure [6] we show the evolution of the field by including the real configuration space 
fields at times t = 0.1, 3.0 and 15.0. The initial field containing intermediate scale structure 
can be seen to develop large scale structures, which unlike ordinary or MHD turbulence 
again appear to be relatively stationary. 

Previous work has considered inverse cascade action in driven, rather than decaying, 2D 
EMHD turbulence js]. There the authors found a forward cascade of energy and an inverse 
cascade of mean squared magnetic potential or anastrophy. Between the forcing lengthscale 
and the artificial dissipative cutoff, the authors found a energy spectrum consistent 

with our results at t = 0.1. It remains unclear though why the spectrum steepens at late 



time. Previous work 



17( 1 has found a tendency toward energy condensation in forced 2D 



classical turbulence. There the condensation is a finite size effect of the biperiodic box 
which occurs after the standard inverse cascade reaches the size of the system. It leads to 
the emergence of a coherent vortex dipole. It is important to note that the dipole contains 
most of the injected energy and since we are simulating decaying EMHD turbulence, we do 
not inject any energy which could power the emergence of such a structure. The real field 
in Figure [6] shows no evidence of coUimated dipole structure and in fact a large number of 
isolated vortices can be seen, a characteristic of fiuid turbulence noted in ^]. Here then 
we deduce we are seeing the first direct demonstration of the dual cascade phenomenon in 
decaying 2D EMHD turbulence. 



IV. CONCLUSIONS 



We have investigated the nature of decaying 2D EMHD turbulence with normal diffu- 
sivity and compared it with classical and MHD turbulence and studies of 2D EMHD with 
hyperdiffusivity. We have found EMHD turbulence experiences an isotropic forward cas- 
cade of energy to higher wavenumber (smaller spatial scales) asymptotically approaching 
increasing Rb (inversely proportional to a dissipation coefficient), in broad 
agreement with previous studies. We have found there is no dissipative cutoff at the pre- 

4/5 

dieted wavenumber k ~ and argue this is consistent with non-local coupling in Fourier 
space, the most important result of this paper. Hyperdiffusivity has previously clouded this 
issue and introduced an artificial cutoff. Only now, when we can avoid its use, has the true 
nature become clear. We have also found that fully developed EMHD turbulence appears 
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to be strongly structured, retaining a similarity to the initial field at late time, very much 
unlike classical or MHD turbulence and a point not noted in previous literature. Our study 
of EMHD turbulence with normal diffusivity has been found to display scale-dependent 
anisotropy in the presence of a uniform background field, in good agreement with previ- 
ous studies employing hyper diffusivity. Further, our results support previous studies which 
found the strength of the anisotropy is linearly related to the external field strength. Finally, 
we have discovered that decaying EMHD turbulence is capable of yielding an inverse cas- 
cade, at least partially transferring magnetic energy from intermediate to large lengthscales. 
This result may be particularly significant for the magnetic fields of neutron stars, where 
the proto-neutron star that emerges from a supernova explosion may well have a primar- 
ily small-scale, disordered field. A Hall-induced inverse cascade may then be a mechanism 
whereby it acquires a large-scale, ordered field. 
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total energy vs. time 




FIG. 3: A plot of energy against time for Rb = 10, 30, 100, 300, 1000 & 3000, increasing from left 
to right. 
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FIG. 4: Fourier power spectra and real space fields at t = 0.2 for 2D EMHD turbulence in the 
presence of a background field. We show 2D Fourier power spectra of a, the poloidal field, for 
Rb = 100 at C = 0, 1, 2 &; 4 across the top row. Across the bottom row we show the corresponding 
Bx fields in real configuration space. 
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Fourier power spectro 
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FIG. 5: Evolution of the Fourier power spectra for Rb = 1000, showing an inverse cascade of 
energy to k < 10. Spectra are shown at t = 0.0 (dashed hne), 0.1, 1.0, 3.0, 6.0 &i 15.0. Note the 
changes in resolution from N = 2048 at i = 0.1 to iV = 1024 at i = 1.0 to = 512 at later times. 
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FIG. 6: Plots of the Rb = 1000 solution in real configuration space at t = 0.1 (top row), t = 3.0 
(middle row) and t = 15.0 (bottom row). The fields have been rescaled onto grids of 100 x 100 
points. 
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